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We propose a new method of calculating nonlinear optical responses of interacting electronic 
systems. In this method, the total Hamiltonian (system + system-pump interaction) is transformed 
into a different form that (apparently) does not have a system-pump interaction. The transformed 
Hamiltonian, which we call the pump built-in Hamiltonian, has parameters that depend on the 
strength of the pump beam. Using the pump built-in Hamiltonian, we can calculate nonlinear 
responses (responses to probe beams as a function of the pump beam) by applying the linear response 
theory. We demonstrate the basic idea of this new method by applying it to a one-dimensional, two- 
band model, in the case the pump excitation is virtual (coherent excitation). We find that the 
exponent of the Fermi edge singularity varies with the pump intensity. 



Nonlinear optics for semiconductors is one of active fields in solid state »hysicsEJ. Recently, research interests of 
nonlinear optics are extended to Mott-Hubbard and charge-transfer insulatorsa, for which electron-electron interactions 
play central roles. The nonlinear optical responses, especially punup-probe spectroscopy, from these systems are 
experimentally studied, and quite interesting phenomena are reportecH. However, it is practically impossible to apply 
the standard method of calculation of x*- 3 -* to such systems, because it is quite difficult to obtain the wavefunctions 
of the ground states and the excited state. 

In this letter, we propose a new approach, in which the total Hamiltonian H (t) (system Hamiltonian H sys + system- 
pump interaction H S y S - pump (t)) is transformed into a different form that (apparently) does not have a system-pump 
interaction. The transformed Hamiltonian is called a "pump built-in Hamiltonian" . If the pump built-in Hamiltonian 
is tractable, one can calculate nonlinear responses (responses to probe_beams as a function of the pump beam) by 
applying the linear response theoryQ to the pump built-in Hamiltonian.El We here present the main idea by explicitly 
constructing the pump built-in Hamiltonian in a two-band model when the pump excitation is virtual (in the sense 
we explain below). It is shown that we can include effects of the pump beam as modification of the parameters of 
the original Hamiltonian in some cases. Namely, the system with a pump beam is equivalent to the system without 
the pump beam with different values of parameters. Using the pump built-in Hamiltonian, we successfully evaluate 
the Fermi edge singularity of the absorption of the probe beam. It is shown that the exponent of the singularity 
varies as a function of the pump intensity. This success indicates that the pump built-in Hamiltonian method could 
have potential to analyze nonlinear optical responses of various systems in which electron-electron interactions play 
important roles. 

We assume that the pump excitation is "virtual", i.e., the quantum- mechanical coherence is preserved when the 
pump beam is shining the sample. The virtual excitation is realized when the spectrum of the pump beam does 
not overlap with the absorption spectrum of the system. If the state vector \^(t)) of the system is initially the 
ground state \G) of H sys before the pump beam is on, it evolves adiabatically the "instantaneous ground state" of 
H(t) = H sys + H sys -p Ump (t) in the rotating frame. Namely, 

|*(t)> = \G[£(t)}), (1) 

where \G[£(t)]) denotes the ground state of H(i) in the rotating frame when t is regarded as a parameter. It is a 
function of the instantaneous value of the envelope function £ (t) of the pump beam at time t. Therefore, it is sufficient 
to calculate the linear response of |G[£]) to a probe beam for each fixed value of £ ; it gives the nonlinear response of 
the system at a time t that satisfies £ — £{t). We denote H (t) for £ = £(t) simply as H. 

To demonstrate the basic idea of the pump built-in Hamiltonian method, we here consider a simple model; a one- 
dimensional, two-band model which has a finite band gap E g . The one-dimensional system may be realized, e.g., 
in a quantum wire made of semiconductors. The two bands, which we call "electron" and "hole" bands, can be 
either the lowest conduction subband and the highest valence subband, or the second and the first subbands of the 
conduction band. The latter case may be more suitable to realize the virtual excitationQ. We assume that the pump 
beam is polarized in the direction normal to the quantum wire direction, whereas the probe beam is polarized parallel 
to the wire. Let Ck,a {c\ a ) and d~k-o {,d_ k _ a ) be the Fermi annihilation (creation) operators for an electron and a 
hole, respectively. We take the model Hamiltonian used in Refs.i and!, which has the following form in the rotating 
frameEHl; 
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H = H 1 + H 2 , 

H r = E (Ec(k)clc ka + E d {k)dl k _ a d. h . a ) (2) 

+ J2 X ( k )( C k<rd- k -<r + ct k _ a cl a ), (3) 



H2 = \H E U i ( l)( C l+ q A'- q a' c k' <y>C ka 



<7,a' k,k' ,q 
~d k q k' 



- 2c l+ q * c k*dl k , +q _ a ,d- h , ) . (4) 



Here, the electron and hole energies in the effective-mass approximation are taken as 

fc 2 

E c {k) = — + E g -u p , (5) 

Em = ^, ( 6) 

respectively, where m c and m d are the electron (c) and hole (d) effective masses, and u> p is the central frequency of 
the pump beam. The \{k) is the product of £ and the interband transition dipole. The U(q) > represents the 
interaction strength between particles. The detuning energy A = E g — uj p is assumed positive and large enough to 
satisfy the condition of the virtual excitation. To simplify the calculations, we assume that the initial state \G) is 
two Fermi seas of electrons and holes with the same Fermi wavenumber k F - This may be realized, e.g., by creating 
electrons and holes by another optical beam; the pump beam is shined after the carriers cool down. To further simplify 
the calculations, we consider the case where 

E g ^> \Ei(k F ) - Ei(0)\ ^> \U\. (7) 

Under these conditions, spontaneous band mixing, such as the instability toward an excitonic insulator, are unfavor- 
able. Moreover, many-body scatterings occur only among electrons and holes near their Fermi points. Furthermore, 
we can neglect interband scatterings (hence such terms have not been included in H). 

Our purpose is to construct Hamiltonian in which the pump term, the second term in Eq. (|^), is eliminated. To 
this end, we first diagonalize Hi by the Bogoliubov transformation 

c klJ = cos 9 k c krT + sm6 k d^_ k _ a , (8) 
d-ka = co<s,6 k d- k -<j - sm9 k cl a , (9) 



as 



where 



H l = Y J E+{k)c\ a c ka + £_(fc)<? fc _ CT <L fe _ CT = Hi, (10) 



E±{k) = ^E+(k) 2 + \{kf ± E-(k), (11) 

cos20 fc = E + (k)/ y / E + (k) 2 + A(fc) 2 , (12) 

sin20 fe - -\{k)/^E+(kf + A(fc) 2 , (13) 

E ± {k) = [E c {k)±E d {k)]/2. (14) 

Note that E + —E c and £L —Ed represent the optical Stark shifts of the (non-interacting) electron and hole, respectively. 
We then rewrite H 2 using the transformed operators, c, c\ d, and d^ . Substituting these operators in H 2 , and ordering 
them in their normal order, we obtain H 2 = H[ + H 2 + H 2 . Here H[ is the one-body part, whereas H 2 and H 2 are 
many-body interactions. The H 2 consists of terms of the forms of cWcc, d^d^dd, and c'd^dc. On the other hand, 
H' 2 consists of other combinations of c, c\ d, and eft . Needless to say, only gauge invariant combinations appear. 
Collecting all the one-body terms, we obtain the one-body part of the transformed Hamiltonian, H\ = H® + H[, in 
the following form: 
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#1 = E (^c(fc)ct£fea + B„(* ; )aL fc _ ff d_ k _ cr ) 

+ ^ A(fe)(c fe(7 <L fe _ CT + dl fc _ CT cL)- (15) 



On the other hand, the form of H 2 is 

cr,a f k,k' ,q 

+^-k-q -J-k'+q -a'd- k > -A-k -a 

-24+g<T £ fc^dLfc/ + g_ CT ^-fc'- CT ')- (16) 

We consider two simplified cases: (a) 17(0) > £/(2fc F )0 and (b) U(Q) ~ U(2k F ). 

Firstly, we consider case (a), for which backward scattering is negligible. In this case, under the condition (0), we can 
take k = ±kp, k' — ±fcp, and q ~ 0. Therefore f/(g, fc, fc') is represented by a single parameter U(0, ±kp, ±fcp) = U ■ 
We then find that H' 2 is negligible. Comparing the explicit forms of H2 and H2, we also find that H2 has exactly the 
same form as H2, and that the parameters are changed as 

E c/d {k) = E +/ _{k) + jsm 2 e k 

-L- S i n 2 kCO82 k + X L S m 2 2e k , (17) 

A(fe) = -^sin20 fe , (18) 

and U — U. Furthermore, the Fermi wavenumber of the transformed system, kp, is equal to fc F because we have 
assumed that i?(fc F ) — E(0)\ 3> \U\ and that the excitation is virtual: in the virtual excitation, the densities of 
"dressed particles" (c and d) are conserved. 

Although our purpose is to eliminate the pump term, a new pump term has been generated through the many-body 
interactions, with the coefficient A(fc). To eliminate this term, we perform Bogoliubov transformations successively: 
since the form of H2 is the same as H2 , we can perform the same procedure successively. The mapping functions from 
{E c , E d , A} — > {E c , Ed, A}, Eqs. ( |l7| ) and Jis|), are denoted by fj(E c , E d , A) (j = 1, 2, 3). Then the recursion relations 
between steps n and n + 1 are written as follows: 

Ei n+ V =f 1 (EW,E d n \\W), (19) 

eW = me<?\eV>,\W), (20) 

A (n+1) _ £g(n) A (n))_ ( 21 ) 

These recursion relations are solved numerically. 
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FIG. 1. The semi-log plot of A' n ^ versus the number of iterations in case (a) !7(0) 3> U(2kp), for several values of 4U/A =6, 
4, 2, and 1 (from the top to the bottom). The A of the original system is fixed. 
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Figure 1 shows the semi-log plot of for several values of U/A, where the pump-system coupling strength of the 
original system, A, is fixed. It is found that A*™-* decreases exponentially with increasing the number of iterations, 
converging to zero for all values of U/A. If we define the decay exponent fi by A( n ) oc exp[— fin], the plot shows that 
the larger the forward scattering strength is, the smaller the decay exponent is. 




#iteration 

FIG. 2. Fhe semi-log plot of A^™- 1 versus the number of iterations in case (a) U(0) S> U(2kp), for several values of 4A/A=6, 
4, 2, and 1 (from the top to the bottom on the left vertical axis). The value 4(7/ A = 1 is fixed. 

Figure2 shows the semi-log plot of A 1 ™) for several values of A/A, where AU/A = 1 is fixed. It shows that A' n J 
converges to zero in the limit of n — > oo for all values of A/A. The plot also shows that the decay exponent is almost 
independent of A of the original system. We thus find that \^ converges to zero after iterations for any values of 
U > and A > 0, and that the decay exponent depends only on the strength of the forward scattering. 

Next the recursion relations for the kinetic energies are considered. Figure 3 shows that the kinetic energies of the 
electron and the hole in the limit of n — > oo, for AU /A — 1 and 4A/A = 1. The curvatures of the dispersions become 
slightly smaller than the original ones (dotted lines). Qualitatively similar results are obtained for other values of 
U /A > 0. We thus conclude that the elimination of the pump term makes the masses of electron and hole heavier. 

Therefore the pump term is eliminated in the limit of n — > oo, and we obtain the pump built-in Hamiltonian in 
case (a). 
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FIG. 3. The dispersion curves of the electron and the hole (black lines) in the limit of n — > oo in case (a) U(0) S> U(2kp)- 
The dotted lines show the initial dispersions. 

Secondly, we consider case (b), where U(0) ~ U(2kp). This situation may be realized, e.g., for a screened interaction 
if the screening length is larger than 1/kp. By dropping the wave number dependence of the interaction strength, we 
find 

cr,a' k,k' ,q 



A 



+d-k-q -J-k'+q -a>d-k> -a'd- k -a 

-2cl+qvCk<Tdl k , +q _ a ,d-k>-<r>y (22) 

and H' 2 = 0. As in case (a), we obtain 

E c/d (k)=E +/ _(k) + ^I 2 

-—I 2 cos26 k + UI 1 $in20 k , (23) 
A(fc) = -LT 2 sin20 fc -Uhcos26k, (24) 

and U — U. The I\t2) are defined by 1\ = (1/2) ^2 k sm29 k and I 2 = Sfc sm2 ^ fc - From successive transformations, 

the recursion relations corresponding to Eqs.([l9|)-(pl|) are obtained. Figure 4 is the semi-log plot of A^ as a function 
of the number of iterations. One can see that A converges to zero exponentially as in case (a), although the decay 
exponent is smaller than that of case (a). The plot corresponding to Fig. 3 is also obtained, which yields the same 
conclusion that the elimination of the pump beam makes the effective masses of electron and hole slightly heavier. 
We thus obtain the pump built-in Hamiltonian in the limit of n — > 00 also in case (b). 
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FIG. 4. The semi-log plot of A'"' against the number of iterations in case (b) U(0) ~ [/(2/cf), for several values of 4(7/ A =6, 
4, 2, and 1 (from the top to the bottom). The A of the original system is fixed. 



Since we have constructed the pump built-in Hamiltonian for both cases (a) and (b), we can discuss nonlinear 
optical responses by applying the linear response theory, the KuboJprmula. For example, the absorption of the probe 
beam in the presence of a pump beam is given by the real part oJi3 



U J -00 



(25) 



where j(q, t) is the current operator, and (• • •} denates the-expectation value in the ground state \G{£]) of the pump 
built-in Hamiltonian. For case (a), following Refs. (B) and (El), we evaluate the power-law singularity of the absorption 
of the probe beam; I {to) ~ \u — Ef] 71 ^^ 1 - It is found that the exponent 77(f) varies as a function of the intensity 
of the pump beam. The explicit form of n{£ ) is obtained by applying the bosonization technique, by which pump 
built-in Hamiltonian is diagonalized in terms of the charge (j = 1) and spin (j = 2) excitations. The result is 



= £ EE 

{e.ti} {1,2} 



(id) 2 



'g)-2u)v) 



(26) 



where [«f,«f|* - TZ{ip 3 )[{v^) 1 ' 2 , {vf) 1 ' 2 ] 1 , = ^(^-)[(«F)" 1/2 , «)~ 1/2 ]', 



5 



c/d _ 



(l/2)vf g] + vpg1. 



± 



(27) 



9i — v f + U/ir, g% = v°p — U/tt, £i = —2U/ir, and £2 = 0. The lZ(<fij) is the 2x2 rotation matrix of angle 
(pj, where 

y c = (1/2) Um-^y/vfvp/ivfgS - vpgf)], (28) 

and if2 = 0. The r){£) is determined by the Fermi velocity and the interaction strength in the pump built-in 
Hamiltonian. In the model under consideration, the Fermi velocity in the pump built-in Hamiltonian, v F , is smaller 
than the one without the pump beam, Vp, whereas the interaction strength is unchanged. Hence, 77(f) > 77, which 
should be observable in pump-probe experiments on quantum wires. Since \E(kp) — E(0)\ 3> \U\ is assumed, the 
difference between v* F and vf is small, hence the change of the exponent is also small. If we take other models, 
however, we expect a larger change. 

Finally, we mention another significance of the pumpj-built-in Hamiltonian method. Sato and coworkers recently 
observed a phase transition induced by a pump beamLLd. Miyashita and coworkers recently suggested that this 
phenomenon can be explained if one assumes that parameters in the Hamiltonian varies as a function of the pump 
intensity!!]. Since the pump built-in Hamiltonian has this property, the pump built-in Hamiltonian method may 
provide for the microscopic foundation of their assumption. However, we do not exclude,other possibilities, such as a 
thermal effect, for the origin of the changes of the parameters in the experiment of RefJi3. Further studies are needed 
for settling it. 
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